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Abstract 

A spherical t-design is a set of points on the sphere that are nodes of a positive equal 
weight quadrature rule having algebraic accuracy t for all spherical polynomials with degrees 
< t. Spherical t-designs have many distinguished properties in approximations on the sphere 
and receive remarkable attention. Although the existence of a spherical t-design is known for 
any t > 0, a spherical design is only known in a set of interval enclosures on the sphere m 
for t < 100. It is unknown how to choose a set of points from the set of interval enclosures 
to obtain a spherical t-design. In this paper we investigate a new concept of point sets on 
the sphere named spherical te-design (0 < e < I), which are nodes of a positive weight 
quadrature rule with algebraic accuracy t. The sum of the weights is equal to the area 
of the sphere and the mean value of the weights is equal to the weight of the quadrature 
rule defined by the spherical t-design. A spherical tg-design is a spherical t-design when 
e = 0, and a spherical t-design is a spherical tg-design for any 0 < e < 1. We show that 
any point set chosen from the set of interval enclosures [15] is a spherical t^-design. We 
then study the worst-case errors of quadrature rules using spherical te-designs in a Sobolev 
space, and investigate a model of polynomial approximation with the Zi-regularization using 
spherical tg-designs. Numerical results illustrate good performance of spherical te-designs 
for numerical integration and function approximation on the sphere. 

Keywords: spherical t-designs, polynomial approximation, interval analysis, numerical 
integration, b-regularization 

AMS suject classifications: 65D30, 41A10, 65G30 


1 Introduction 

For a d-dimensional sphere 

S'i := {x = (xi,... ,Xd+i)'^ G I ||x||^ = 1}, 

where || • ||2 means the Euclidean norm, a spherical t-design m for a given positive integer t is 
a set of N points Xn = {xi,... ,XAr} C such that 

N 
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holds for all spherical polynomials p with degree < t, in which |S'^| is the area of 8 “^ and dw^ 
is the surface measure. For a numerical integration rule on the sphere, we say that the rule 
has algebraic accuracy t, or the rule is of t-algebraic accuracy, if the rule integrates all spherical 
polynomials with degree < t exactly. A spherical t-design establishes a positive equal weight 
quadrature rule with algebraic accuracy t for numerical integration on the whole sphere, which is 
also proved to perform well for numerical integration of spherical functions belonging to Sobolev 
spaces, see [ 12 ] for detail. 

The existence of spherical t-designs for arbitrary degree t was proved by Seymour et al 
|26| in 1984. Consequently, a natural problem is to find the minimal number of points such 
that ([H) holds for fixed t and d, denoted as N{d,t). In 1993, Korevaar et al [22] proved 
that N{d,t) < and conjectured that N{d,t) < Cdt'^ for a sufficiently large positive 

constant Cd depending only on d. This conjecture was then proved by Bondarenko et al [TO] 
in 2011. For d = 2, there is an even stronger conjecture by Hardin et al m saying that 
N{2,t) < + as t — oo. Numerical evidence supporting the conjecture was also given in 

[IS1E7]. Spherical t-designs have been extensively studied from various viewpoints, among which 
the application to polynomial approximation and the number of points needed to construct a 
spherical t-design have been paid great attention, see [siiaisiEiiTiisiiiaiiniisaisTj. Moreover, 
numerical methods have been developed for finding spherical t-designs. Usually in those methods 
finding spherical designs is reformulated to problems of solving a system of nonlinear equations 
or optimization problems, see [ailHlET]. Numerical results suggest that those methods can hnd 
approximate spherical t-designs with high precision. 

However, in general, for any given t the exact location of a spherical t-design = 
{xi,..., xat} is unknown. The best we know is that there is a set of points x*, i = 1, ..., N 
such that a set of narrow intervals defined by 

Xn = {[x]i = C(xj,7i), i = Xj G 7i > 0} C 

can be computed to contain a spherical t-design for d = 2, iV = (t + 1)^ and t < 100 in [15], 
where 

= {x G I cos“^(x • Xj) < 7 ^}. 

Among all the spherical polynomials with degree < t, if the zero polynomial is the only 
one that vanishes at each point in the set Ajv G S^, then the point set G is said to be 
fundamental with order t. In 2011, Chen et al [15] proposed a computational-assisted proof for 
the existence of spherical t-designs on with N = {t + 1)^ for all values of t < 100. An interval 
arithmetic based algorithm is proposed to compute a series of sets of polar coordinates type 
interval enclosures containing fundamental spherical t-designs. By choosing the center points 
of each interval enclosures, an approximate spherical t-design can be obtained and numerical 
results show that the Weyl sums of these point sets are very close to 0. However, though we 
have known the existence of spherical t-designs for t < 100 and N = {t + 1)^ in a set of interval 
enclosures, we still can not obtain an exact spherical t-design in the set of interval enclosures for 
a positive equal weight quadrature rule with algebraic accuracy t. Motivated by this problem, 
in this paper we relax the equal weights to the ones whose mean value is still but can be 
chosen in an interval with respect to a number 0 < e < 1 . 

Definition 1. (Spherical te-deisgn) A spherical te-design with 0 < e < 1 on is a set of 
points A^ := {xj,...,x^} C such that the quadrature rule with weights w = {wi, ... ,wn)'^ 
satisfying 

i = l,...,A, 


( 2 ) 
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is exact for all spherical polynomials p of degree at most t, that is, 


N 

V Wip{xi) = / p(x)dwrf(x) 


(3) 


The concept of spherical te-designs establishes a bridge of spherical t-designs and ordinary 
positive weight quadrature rules with algebraic accuracy t. 

Remark 2. A spherical t-design is a spherical tQ-design with e = 0. By letting p{x) = 1 in (0) 
we can obtain Wi = |S'^| and thus 0 < Wi < |§'^| for i = 1, ... , N . 

Since the existence of spherical t-designs has been proved for arbitrary t, and a spherical 
t-design is also a spherical t^-design for arbitrary 0 < e < 1, we have the existence of spherical 
t^-designs. Due to the relaxation of the weights w, an important advantage is that we can have a 
positive weight quadrature rule with algebraic accuracy t. Moreover, our numerical experiments 
show that with the increase of e we can get numerical integration with polynomial precision 
using fewer points than spherical t-designs. 

The rest of this paper is organized as the following. In Section 2, we discuss the relationship 
between spherical t^-designs and spherical t-designs when they are fundamental systems with 
same number of points. Based on the results we study the sets of interval enclosures containing 
fundamental spherical t-designs in Section 3. We prove that all point sets arbitrarily chosen 
in those sets of interval enclosures computed in m are spherical t^-designs. In Section 4, 
we analyze the worst-case errors of spherical t^-designs for numerical integration on the unit 
sphere S^. Numerical results show that compared with the equal weight quadrature rules using 
spherical t-designs, the worst-case errors can be improved with relaxing the weights to define 
quadrature rules using spherical te-designs. In Section 5, we investigate an I 2 — h regularized 
weighted least squares model for polynomial approximation on the two-sphere using spherical 
t^-designs and present numerical results to demonstrate the efficiency of the I 2 — h model. 

In this paper we concentrate on the case d = 2. Throughout the paper we assume that all 
the points in a point set on the unit sphere are distinct. The computation is implemented in 
Matlab 2012b and done on a Lenovo Thinkcenter PC equipped with Intel Core 17-3770 3.4G Hz 
CPU, 8 GB RAM running Windows 7. 


2 Spherical ig-designs: neighborhood of spherical ^-designs 

In this section we will study the relationship between spherical tg-designs and spherical t-designs 
when they are both fundamental systems and have the same number of points. A spherical t- 
design defines an equal weight quadrature rule with algebraic accuracy t while a spherical 
te-design defines a positive weight quadrature rule with algebraic accuracy t. Based on these 
two properties, in this section we study a neighborhood of spherical t-designs. Denote 

Pt := Pt(§^) = {spherical polynomials of degree <t} 

= spanjy^^fc : i = 0,... ,t, k = 1,... ,2i + 1}, 

as the space of all spherical polynomials with degree < t on the unit two-sphere Here 
is a fixed L 2 -orthogonormal real spherical harmonic of degree i and order k, which means 


[ 


Yi^k{x)Yi:^k:{x)duj{x) = 5i^ii5k,k', = 0, ■■■ ,t-, k,k' = l,...,2i + l, 


(4) 
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where da;(x) = da; 2 (x), and is the Kronecker delta. It is well known that 

t 

dt :=dim(Pi) =^(2£ + l) = (t + l)2. (5) 

i=o 

By the addition theorem we have 


2£+l 


ye,k{^)ye,k{y) = • y)> 


k=l 


47r 


( 6 ) 


which implies J2k^=i ylki^) = where Pi, £ > 0 denotes the Legendre Polynomial and x • y 
denotes the Euclidean inner product. Hence we obtain 


/ 21+1 

ll^£,fcllc(§2) < max f Y ylki 

\ k=l 


X 


2i+l 

Att 


loi k = 1,... ,2£ -\-1, £ > 0. 


(7) 


Indeed, a spherical coordinate form of real spherical harmonics can be represented as (see mm) 

V2Ni^kPi^^~^ {cos 9) cos k^p, k = 

Ni^kPe {cos 9), k = £ + l, (8) 

y/2Ni^j.PY^~^ {cos 9) sin k(p, k = £ + 2,... ,2£ + 1, 


Y£,fe(x) = < 


with 


X = 


sin 9 cos ip 
sin 9 sin (p 
cos 9 


where O<0<7r, 0<(/9<27r and Ni^k are the normalization coefficients 


Ni^k 


2£ + l{£-\k-£-l\)l 
47r {£ + \k-£-l\)L 


k = l,...,2£ + l. 


When taking k = £ + I and 0 = 0 we can obtain Y£^£+i(x) = y with x = (0,0,1)"^. 
Therefore, d?]) is a sharp upper bound of ||k£,fc||c(§ 2 )- With the fact that all spherical harmonics 
are a basis of P* we have that a finite point set Xjv = {xi,..., xat} is a spherical t-design if and 
only if the Weyl sums 


N 

^F^,fc(x,) = 0, k = l,...,2£ + l, £ = !,...,t, (9) 

hold, see [IZIET] for details. Let Xj\f = {xi, • • • ,xjv} C and X'j^ = {x'^, • • • ,x(y} C be 
two point sets on the sphere. To describe the relationship among points and point sets, we 
introduce the following definitions of distances. 

1. Define the geodesic distance between two points :x-i,Xj G as 

dist(xj,Xj) = cos“^(xj • Xj). 
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2. Define the separation distance of a point set as 

p{Xm) = mincos“^(xj • Xj), 

which represents the minimal geodesic distance between two different points in 

3. Define the least distance from a point x G to a point set X^ as 

dist(x, Xjv) = niin dist(x, Xj) = min cos“^(x • x*), 
i i 

which can be seen as the geodesic distance from x to its projection on X^r. 

4. Define the Hausdorff distance between two point sets and X^ by 

a{Xj\f, X'^) = max{maxdist(x', Xtv), maxdist(xj, X^)} 

i i 

= max{maxmincos~^(x'• Xj),maxmincos~^(x^ • Xj)}. (10) 

i j j i 

Note that a{XN,X'^) = a{X'j^, Xn) and a{XN, X'j^) = 0 if and only if Xn = X'j^. 

Remark 3. For two point sets X^r and X'^, if a{Xi^,X'^) < ip(Xjv), then for each Xj G Xat 
there exists a unique x' G X'^ nC(xj, where 

C{xi, ^p{Xn)) = {x G I cos"^(x • Xj) < ^p(XAr)}. 

On aceount of the relationship between C(xj, ^^(Xtv)) and X'^, in what follows we will de¬ 
note x( as the point belonging to C(xi, ^p(X]y)) for sake of consistency under the assumption 
a{XM,X'j,) < \p{Xm). 

For a point set Xjv = {xi, • • • ,X 7 v} C we define the matrix Y{Xjq) G with its 

elements as 

Yj_£2+fe(XAr) = Y);^fc(xi), i = 1,... ,X, /c = 1,... ,2^ +1, £ = 0,... ,t. (11) 

Note that |§^| = dvr. An equivalent condition of spherical t^-designs given in [T3] states that a 
point set X^ := {x},..., x)y-} C is a spherical t^-design if and only if 

Y(X^)^w - = 0 and < w < ^ e, (12) 

where ei = (1,0,... , 0)^ G and e = (1,..., 1)”^ G 

For two matrices constructed by two near enough point sets on we have the following 
property. 

Proposition 4. For any two point sets X^, X)^ C satisfying a{XN,X'j^) < ^p{Xn) there 
always holds 

||(Y(X^)-Y(X)v))^||oo = ||Y(X^)-Y(X)v)||i 

< N{t + 1)^ a{XN,X'j^), (13) 

where Y(XAr), Y(X)y) G are matrices defined by 071) which depend on t. 
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Proof. For a point x* G Xn, by Remark[3]we let x' be the unique point located in C(xj, ^p{Xn)) 0 
X'j^. Let be the restriction of on the great circle through these two points. Then ^ 
is a trigonometric polynomial on the sphere and by Bernstein’s inequality m and dZD we obtain 

|h€,fc(xj) - h£,fc(x-)| = \Qe,k{^i) - Qe,k{^'i)\ 

< COS“^(Xj • x') sup \ 

< COS“^(Xi •x')(t + l)sup|(5£,fe| 

< cos"^(xi • x-)(t + l)||Y£,fc||c(S2) 

/2/ -I- 1 

< a{X^,X'^){t + l)^^, (14) 

where the last inequality is obtained by ([71). Together with (|13l) and da we have 

N 

||Y(X^)-Y(X)v)||i = max V |Y,,fc(x,) - Y£,A:(x')| 

J=1 

< N{t + l)\J a{XN,X'^). 


□ 


Let X^ = {x)*,..., x^} C be a fundamental spherical t-design. Given a number a* > 0, 
denote the neighborhood of Y)y with radius a* by 

C{X%,a*) = {Xn C §2 : < a*}. 

The following lemma indicates that any point set contained in a small enough neighborhood of 
a fundamental spherical t-design is a fundamental spherical t^-design. 

Lemma 5. Let be a fundamental spherical t-design with order t and N = {t + 1)^. Then 
any point set X]\j G C{X^,a*) is a fundamental spherical t^-design with 


Ta*\\Y{X%)-^\ 


1-tu*||Y(X^ 


i-ill. - 


< e < 1, 


where 


a < — mm 
2 






NJ 5 


(15) 

(16) 


x\\X{X% 

with T = yj‘^{t + lf. 

Proof. This lemma can be proved by showing that for any point set Xjy £ C{X^, a*) we have 

Y(X7v)^w = (17) 


in which Y(X 7 v) is nonsingular and ||w — 


< 


477 

N • 


From Proposition m for any point set Xn with cj(XAr,X^) < a* we have 
||Y(X^)-i||i||Y(X^)-Y(X^)||i < ||Y(X^)-i||ir^7* < 1. 


Hence Y{X]\f) is nonsingular, that is, Xjv is a fundamental system with order t. 
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By the fact that is a fundamental spherical t-design and Q we have 

^Y{X%fe = V^e,. (18) 

Together with the well known perturbation theorem of linear systems in m Theorem 2.3.9, 
pp.l35] and 

\\l-Y{X<^))-^Y{XMf\\^ < ||Y(X^)-i||i||Y(X^)-Y(X^)||i < 1, 


we have 


-felloo ^ ||(Y(Y^)^)-i||oo ||(Y(Yjv)-Y(X^))^|U 

^ello, - l-||(Y(Y^)^)-i||oo||(Y(Y^)-Y(Y^))^||oo 


Therefore, we obtain 


47r 11 47r ||w 

" 'N ll^elU 

< dTT ||Y(Y^)-ii||Y(Y^)-Y(XA,)||i 

- Y 1 - ||Y(Y^) - Y(Ya,)||i||Y(Y^)-i||i 

< 47r ro-*||Y(Y^)-^||i ^ 

- Yl-rfT*||Y(X^)-i||i Y‘ 


(19) 


Hence the vector w is positive which implies that Yjv is spherical te-design with e satisfying 

([IS]). □ 

Based on above lemma we can deduce the following result which describes an upper bound 
of the radius of the neighborhood of a fundamental spherical t-design, in which any point set 
located in the neighborhood is a fundamental spherical 4-design. 

Corollary 6. Let X^ be a fundamental spherical t-design with order t and N = [t + 1)^. For 
any 0 < e < 1, if 

then Xpi is a fundamental spherical t^-design. 

Proof. By the fact that a{XN,X^) < a* which is defined in (USD, we have Y (X^) is nonsingular 
and 

II _4^ II ^47r ra(Yjv,X^)||Y(Y^)-i||i dvr 

||w ^e||oo_ 

Hence, from (12011 . we derive 


, .dvr , .dvr d7r(l — e) ^ 

(1 - e)-e < w < (1 + 6)-e < ^ ^ ^ e. 


We complete the proof. 


□ 
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3 Interval analysis of spherical t-designs and spherical t^-designs 

In this section we will study the sets of interval enclosures containing fundamental spherical 
t-designs. In the last section we describe a neighborhood of a fundamental spherical t-design in 
which any point set is a fundamental spherical tg-design. A series of sets of interval enclosures 
of spherical t-designs can be computed in [15] , but the exact location of spherical t-designs can 
not be obtained. In the following we will show that any point set in the set of the interval 
enclosures given in m is a fundamental spherical tg-design. Let 

Xat = {[x]j = C(xj, 7 j) C i = l,...,N} (21) 

be a set of spherical caps, with Xj as the center point and 7 * as the radius. Additionally, let 
Xn = {xi,... ,X 7 v} C be the set of center points. Define the radius of Xjv by 


rad(Xjv) = max 7 *, 

l<j<7V 

and the separation distance of Xat by 


p{X]\f) = min dist(xj,Xj). 

^ # j 

Xj e [x]j,Xj G [x]j, 


We say that X^r is an interval enclosure of a point set A^r = {xi,..., xat}, denoted as A^r G X^r, 
if Xj G [x]i and Xj ^ [x]j for i = 1,..., A and i / j. 

Assumption 7. LetXjv defined by i21\) be a set of intervals. Assume that 

1. there exists a spherical t-design G Xjv; 

2. Y(Ajv) is nonsingular. 

The assumption that Y(Ajv) is nonsingular also implies that Y(Ajv) is a square matrix 
which requires A = (t + 1)^. In the following theorem we show that under Assumption [7] if 
rad(Xjv) is smaller than a certain number, then Y(Ajv) is nonsingular and (I20p holds for any 
Ajv G Xjv. 


Theorem 8 . Under Assumption^ any point set X^ G Xjv is a fundamental spherical t^-design 
with 

^ \-1 I , 

1 


2Trad(Xjv)||Y(Ajv)-^' 


1 -4Trad(Xjv)||Y(Ajv) 


-i| 


< e < 1, 


if 


rad(Xjv) < min ( \p(X]\f), - - -^- 

U ^ 2(l + 2e)r||Y(Ajv)-i||i 


( 22 ) 


( 23 ) 


Proof. For any Ajv G Xjv it can be concluded that 


p{Xn) > pO^n)- 

Hence for any two point sets Ajv, A)y- G Xjv we have 

(t(Ajv,A^)< max max dist(xj, yj) = 2rad(Xjv). 

l<i<N Xi,yie[x]i 
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First we show that Y{X]\f) is nonsingular for any G Xtv- By Proposition 0] we have 

||I-Y(X^)-^Y(XAr)||i < ||Y(X^)-i||i||Y(Y^)-Y(Y,v)||i 

< ||Y(Y^)-i||ira(Y^,XAr) 

< 2r||Y(X7v)"iirad(XM) < 1. 

Hence all the point sets X]y G X^r including X^ are fundamental systems with order t. More- 


aiXN,X%) < -p{X%). 


over, from (|23]l we have 


By ([23]) we can also have 

2(1 + e)rrad(Xjv)||Y(YAr)-^||i < e - 2 errad(Xjv)||Y(Y 7 v)-i||i. 

By Corollary 2.7 in [29l pp.ll9] it can be concluded that for arbitrary Xjq G Xtv we have 


|Y(Y^) 


-ii 


Ii < 


|Y(Y^)-i| 


1 - ||Y(Y^)-i||i||Y(Y^) - Y(X,v)||i 
Then together with Proposition 0] we have 

a{XN,X%) < 2rad(X^) 


(24) 


< 

< 

< 


e(l -2Trad(X7v)||Y(Yjv) 


-i| 


1 


(l + e)r||Y(Y^)-i||i 

6(l-||Y(Yjv)-i||i||Y(YA,)-Y(Y^)||i) 

(l+e)r||Y(YAr)-i||i 

e 


(i+6)T||Y(x^)-i|ir 

By Corollary [6] we have that any point set X^ G X^r is a fundamental spherical tg-design. 
Additionally, together with 


dvr 

|w-—ello, < 


< 


dTT ||Y(Y»,)-^||i||Y(Y)^)-Y(Xa,)||i 
Yl-||Y(Y^)-Y(X^)||i||Y(Y^)-i||i 
dvr 2rrad(XAr)||Y(X^ 


0 ^-l| 


and 


we have 


l-2rrad(X^)||Y(Y^)-i| 
2rrad(X^)||Y(Y^)-ii 


l-2rrad(X7v)||Y(X^) 
2Trad(X7v) 


-i| 


l|Y(Y^)-^| 


< 


< 


1 - ||Y(Y^)-i||i||Y(Yjv) - Y(YAr)||i 


1 — 2Trad(X7v) 


l|Y(Y^) 


-ii 


1- ||Y(XAr)-i||i||Y(XAr)-Y(YAr)||i 


-i| 


2Trad(Xjv)||Y(Yjv) 
l-4Trad(XAr)||Y(YAr) 


-i| 


Then we have that any point set G X^r is a spherical tg-design with 

2rrad(X^)||Y(Yjv)-ir 


1 -4Trad(XAr)||Y(Y7v) 


-i| 


< e < 1, 


under the assumption of rad(XAr). 


□ 
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Theorem [8] proves that an arbitrarily chosen point set in a set of interval enclosures of 
a fundamental spherical t-design is a spherical t^-design if rad(XAr) is smaller than a certain 
number. In the theorem we discuss the interval enclosures defined by spherical caps as [x]j = 
C(xj, 7j) = {x G §^| cos“^(x-Xj) < 7j} for which many nice properties of real spherical harmonics 
can be adopted. However, in practice, to reduce the spherical constraint of points and the 
dimension of variables, the spherical coordinate form of the points are preferable to compute 
the interval enclosures, see [15]. For a point Xj G C S^, denote 6i, cpi as its spherical 
coordinate. Then in m a series of intervals [0]i = [0i,9i], [p]i = [(pi,ipi] are computed such 
that Zat = {[z]i,..., [zJtv} is a set of interval enclosures of a well-conditioned spherical t-design 
[2], in which each element in Z^r is defined by 


z i = 


/ sin([6']i)cos([¥5]j) ^ 
sin([6»]i)sin([(^]i) 

V cos([6i]i) / 


i = l,...,N. 


(26) 


In this sense, different from the interval enclosures defined by the spherical caps, each interval 
enclosure computed in |15] is a rectangle as [0]i x ['p\i. Therefore, there remains a gap between 
real computation of interval enclosures of spherical t-designs and our analysis above. Naturally, 
a strategy to overcome this gap is that for each spherical rectangle in |15] we construct a 
spherical cap which is as small as possible to cover it. For the spherical rectangle \9]i x [{p\i = 
[9i-,9i] X [ipi,(pi] its four vertices can be written as 


/ sin(0i)cos(^i) \ 


/ sin(0j)cos((^j) \ 

sm{9i)sm{(pi) 

, = 

sm{9i)sm{pi) 

0 

0 

C/3 


V cos{9i) ) 

/ sin(0i) cos((^i) \ 


/ sm{9i)cos{(pi) \ 

sm{9i) sm{ipi) 

; ^2,4 — 

sm{9i)sm{pi) 

V cos{9i) j 


V cos((9i) ) 


It can be shown that there exists a point Xj defined by [9, ip\ G [9]i x [(/?]j satisfying 


dist(xi, Xjj) = dist(xj, Xj^fc) for j, A: = 1, 2,3,4 


(27) 


and [z]j C C(xj, 7j) with 7* = dist(xj, Xj^). However, computing such point Xj is time-consuming 
and imports large round-off errors when the radii of interval enclosures are small. Instead of 
computing Xj, we investigate another strategy to compute the spherical caps to cover spherical 
rectangles which is coarser but more practical. For a spherical coordinate interval [z]j, we use 
the center point of the interval [0j, 9j\ x \pi, (pi] to dehne a point as 


/ sin(i(6'j0i))cos(i((/?j-h ^i)) \ 


Xj 


sin(i(6'j -h 9i)) sin(i(y5i -h pi)) 


(28) 


V 


cos(i(6'i 9i)) 


/ 
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Note that the spherical coordinate of Xj is the center point of the interval [9]i x \Lp\i but itself 
is not necessary to be the center point of [z]j in the form of the spherical coordinate. Still, we 
have that 

dist(xj,Xj,i) = dist(xi,Xi,2), dist(xj, Xj^a) = dist(xj, Xj,4). (29) 

It is obvious to obtain that the distance between Xj and any point in [z]j does not exceed the 
maximum of the four distances in (1291) . Therefore, if we let 

7j = max{dist(xj,xj,i) , dist(xj,Xj^a)}, (30) 

then we have 

[z]i C C(xi,7j). (31) 

Consequently, we call the set of spherical caps 

Xat = {C(xj,7j)} 


as a cap-cover o/Z^v with x*, 7* defined in (1251) and (ISUIl . Similar with set of spherical caps Xtv, 
we define the radius and separation distance of Z^v by 

rad(ZAr) = max { max{dist(xj, Xj^) , dist(xj, Xj^a)}} = rad(XAr), (32) 


and 

p(Zn) = min {dist(xj, x^) - 7* - 7^} = p(XAr). (33) 

« 7^ 3 

1 < i,j < N 

And then we have the following corollary for the lower bound of e for Z^v- 


Corollary 9. Let "Ln be a set of spherical rectangle and its cap-cover X^v satisfy Assumption 
Then any point set Xn G Ztv is a fundamental spherical t^-design with 


2rrad(Xjv)||Y(Xjv)-^||i 

l-4Trad(Xjv)||Y(A^)-i||i 


(34) 


if 

rad(Z7v) < min ( -pCLj^), -^—=- ) . (35) 

V4 ^ 2(l + 2e)T||Y(A^)-i||iy 

The proof of the corollary follows the same manner of Theorem [5J 

Based on Corollary O the lower bounds of e, denoted by e, for the sets of interval enclosures 
provided in m can be computed. The data containing the sets of interval enclosures for the 
parameterization of the spherical t-designs and relative programs can be downloaded from the 
website 

http://www-ai.math.uni-wuppertal.de/SciComp/SphericalTDesigns. The computational results 
are shown in Fig 3.1 and Table 3.1. 

In Fig. Ol, we report the lower bounds of e, denoted as e, for sets of interval enclosures 
computed in [15] for t = 2,, 100 (for t = 1 we have known that the regular tetrahedron is 
a spherical t-design so that e = 0 and we would not consider this case here), based on formula 
(1341) . In this hgure we also plot a function 




(36) 


to approximately estimate the track of e with respect to t. From the figure we can conclude 
that the lower bound of e grows with the increase of t in an order about 6.9. Additionally, by 
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Figure 1: e for i = 2,..., 100 

(1341) it is known that the norm of Y(Xjv)~^ is also very important in the process of estimating 
e. Fortunately, since the sets of interval enclosures computed in m seek to include well- 
conditioned spherical t-designs [2], the growth of lower bounds of e keeps stable for all the t 
considered here. 


Table 1: Information for sets of interval enclosures Zjv for some selected t 


t 

rad(ZAr) 

p{Zn) 

e 

e for Xn 

10 

1.843454e-12 

3.396362e-01 

6.748890e-08 

6.694645e-14 

20 

1.515848e-ll 

1.805783e-01 

5.480524e-06 

1.783018e-13 

30 

5.588085e-ll 

1.249714e-01 

7.888659e-05 

2.480238e-13 

40 

1.044163e-10 

9.203055e-02 

4.043164e-04 

5.339063e-13 

50 

2.199182e-10 

7.638945e-02 

1.862348e-03 

5.057066e-13 

60 

4.006638e-10 

6.302748e-02 

6.502352e-03 

6.747935e-13 

70 

6.143914e-10 

5.421869e-02 

1.820130e-02 

8.820722e-13 

80 

1.220430e-09 

4.771142e-02 

6.050880e-02 

1.151368e-12 

90 

2.089473e-09 

4.264961e-02 

2.066649e-01 

1.228462e-12 

100 

2.273791e-09 

3.846343e-02 

4.420562e-01 

1.880540e-12 


We also report some information of the interval enclosures and their theoretical lower bounds 
of e for some selected t in Table 3.1. Not only e but also the radii and separation distances of 
are also shown in the table. We can see that the radius of each interval enclosure is far smaller 
than their separation distance, which means that the assumptions in the above lemmas and 
theorems are satisfied. For a fixed t, the set of center points for each interval ([0]i, [<p]i) £ Zjv, 
i = 1, ... ,N, denoted as Xn = {xi ... ,xn}, with x* defined by (I25|) . can be regarded as an 
approximation of a fundamental spherical t-design. For each new point set Xjsf we compute e 
by (1341) so that the equality holds exactly. As shown in the table, the values of e for each Xw 
are all very small positive numbers, and grow with the increase of t. This means that Xw which 
is selected properly from the interval enclosures is a spherical fg-design. 
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Remark 10. To obtain an approximate spherical t-design as accurate as possible, in [15] the 
radius of the set of interval enclosures rad(Z 7 v) is computed to a small scale around 10“^*^. With 
the introduction of the concept spherical t^-designs, it has been shown that any point set selected 
in Ztv is a fundamental spherical t^-design with small enough rad(ZAr) and e is increasing in 
rad(Z 7 v). As a result, to reduce the difficulty of computing Z^v, one may relax rad(ZAr) to a 
larger scale with keeping e < 1 holding. Then all point sets selected inZ]\f are still fundamental 
spherical t^-designs with 0 < e < 1. 


4 Worst-case errors of quadrature rules using spherical tg-designs 


In this section we will investigate the worst-case errors for quadrature rules with algebraic 
accuracy t using spherical t^-designs in Sobolev spaces which are finite-dimensional rotationally 
invariant subspaces of (7(8^). The bizonal reproducing kernel will be used in the analysis, which 
has been wildely applied to analyze approximations on the sphere [T2l [T3l [201 EH 131] . About 
equal weight quadrature rules, Brauchart et al [T2| recently develop a way to compute their 
worst-case errors in Sobolev spaces. In this section, we intend to extend their method to non¬ 
equal but still positive weight quadrature rules and show the performance of spherical tg-designs 
in numerical integration. 

In this section we follow notations and definitions from [12] . Denote the space of square 
integrable functions on by L 2 (S^). Then it is a Hilbert space with the inner product 

(/> 9)h2m = ^2 f,g e L 2 (§^), (37) 

and the induced norm as 

II/IIl2(s2) = (^^J/(x)|^dw(x)^ , /gL2(§^). (38) 

The Sobolev space BI^(§^) can be defined for s > 0 as the set of all functions / G L 2 (§^) with 
whose Laplace-Fourier coefficients 


fe,k = (/,T£,fc)L .§ 2 ) = / /(x)Y)-,fc(x)da;(x), 

satisfying 

oo 2i+l 2 

X] hk < 

e=o k=i 

where + 1). Obviously, by letting s = 0 we can obtain ]H‘^(§^) 

norm of ]HI*(§^) can be defined as 


oo 2£+l 


~flk 


]_l=Q k=l 

where the sequence of positive parameters af'^ satisfies 

~ (1 + A,)-* ~ (^ + l)-"7 

Correspondingly, the inner product of BI^(§^) can be defined as 

oo 2£-|-l ^ 

(/)5)h'> = X/ X/ ~(fjfhk9e,k- 

£=0 k=l 


(39) 

(40) 

= L 2 (§^). Then the 

(41) 

(42) 


( 43 ) 
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For a point set Xjy and a weight vector w, we define the numerical quadrature rule and the 
integral of a function / on as 


^ r 

Q[XN,w]{f) := ■= y^^/Wdw(x), 


(44) 


The worst-case error of the quadrature rule Q[X7v, w] on EI'^(§^) can be defined as mm 


Es{Q[Xn,^]) :=sup{|Q[XAr,w](/)-/(/)| :/€]H*(§2 ),||/||h= <1}. (45) 


The Riesz representation theorem and the additional theorem assure the existence of a repro¬ 
ducing kernel of 


Ks{^,y) = '^{21 + l)af^Pei:xi ■ y) 
e=o 

oo 2£+l 

= 4*^^bfcW^bfc(y)- (46) 

£=0 k=l 


Together with the property of reproducing kernel Ks{-,-) defined in ()46p and the addition 
theorem, it is shown in m that 


(i?,(Q[XAr,w]))' = 


n 2 


sup \Q[XN,w]{f) - I{f)\ 


f 6 

II/IIh'! < 1 
N 


^ |^^i^s(-,x) - / R:s(-,x)da;(x) 
i=i 


Assume that Q[A7 v,w] has algebraic accuracy t. Then with the following equality 



Ks{-k, •)da;(x) 


= a, 


(^) 


the worst-case error could be reformulated as 


{E,{Q\X^M)f 


oo 2£+l 


E E 


£=l k=l 

\j=i / 


N N 


EE 


WiWj 

IGvr^ 


oo 2£+\ 


£=1 k=l 


(47) 


Reproducing kernels for ]HI®(S^) for s > 1 can be constructed utilizing powers of distances, 
provided the power 2s — 2 is not an even integer. Indeed, it is known (cf., e.g., OHS]) that the 
signed power of the distance, with sign (—l)^"*"^ with L := L{s) := [s — Ij, has the following 
Laplace-Fourier expansion 


(_l)i+i|x - y|2--2 = (-l)^+iF 2 _ 2 .(§') + ^aj*^(2£ + l)P,(x • y), (48) 

£=l 
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where 


and 


Thus we have 


^2-2.(§') := 



|X - y 


2s-2 


da;(x)da;(y) = 2 


r(3/2)r(s) 


0rr(l + s)’ 


As). y p 

-^2-2.(S)-^ 


> 1, 


(1 ~ s)e _ r(i + fi) T{i +1 - s) ^ r(i - s) „_ 2 s ^ p- 2 s 
(1 + s)j> r(i — s) r(i’ +1 + s) r(i + s) 


(49) 

(50) 


(-l)^+2(^^_^^(g2)_|^_y|2.-2) ^ ^a{^)(2£ + l)P,(x-y) 


i=i 

oo 2^+1 

H) 


^^T(x)^i’,fc(y)- 


e=i k=i 


(s) 

Note that for ' we have 


r(i)r(s) 


(s) ^ r,2s-l_ 

^ 0F(-i)^+ir(i + s) 


£ as £ ^ oo, 


(51) 

(52) 


(53) 


and when 1 < s < 2, which means L = L{s) = 0, we have > 0 for all £ = 1,.... Therefore, 
we regard the left hand side of (1511) as the reproducing kernel of ]HI*(S^), which is 


P:^(x,y) = V 2 - 2 s(§^) - |x - yp^ 


and then we obtain 


N N 

{Es{Q[XN,w])f = EES (142-2s(S^) - |xi - Xjl^'' 2). (54) 

i=i j=i ^ 

(5) 

For the case s > 2, we know that a) ' > 0 does not hold for all £ = 1,.... In this situation, we 
let 

iF,(x,y) = (1 - (-l)^+i)y 2 - 2 s(§') + Ql(x • y) + (-l)"'+'|x - yp-^ 

with 

L 

Qi(x • y) := ^((-1)^+1-^ - l)a^;\2e + l)P,(x • y), x,y G 
i=i 

which changes the signs if the negative coefficients in (1481) . Hence the worst-case error on 
IHI^(§^) with s > 2 can be represented as 


N N 

{Es{Q[XN,w])f = J2Y1 jS {QU^i • Xj) + (-1)^+Vi - Xj 

i=l j=l ^ 

-(-1)^+V2_2.(S2)) . 


2s-2 


(55) 


In what follows we will compute the worst-case errors of quadrature rules using spherical t^- 
designs with algebraic accuracy t. In this experiment we choose e = 0.1 for spherical tg-designs 
and use (1121) to find a spherical p-design, which is a system of nonlinear equations. The system 
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can be solved by minimizing its least squares form using a smoothing trust-region filter method 
proposed in [T3]. Note that the number of points needed for constructing spherical te-designs 
may decrease with the increase of e. Thus in the computation of spherical tg-designs we always 
attempt to find the one with a possible minimal number of points, denoted as N{t,e). The 
detailed process for finding spherical 4-designs can be found in [2]. In the numerical test of 
computation of spherical fg-designs it is found that a possible minimal number of points satisfies 
[(t-|-1)^/3] -|-1 < N{t,e) < [(t-|-2)^/2] -|-1. In the numerical test in current and next sections, 
the spherical tg-designs are chosen with N = N{t,e). 


[ X spherical t -design | 



O spherica] t-design 

O spherical t-design ] 

10 ’ 

X O 

X spherical t -design 


V 


'X 


(a) s = 1.5 


(b) s = 5.5 


Figure 2: Worst-case errors for spherical to.i-designs and spherical t-designs 

The worst-case errors of quadrature rules using spherical to.i-designs in ]HI®(S^) for s = 1.5 
are illustrated in Fig. 2(a), For comparison, the worst-case errors for quadrature rules using 


approximate spherical t-designs computed in [27] will also be implemented. For all spherical 4.1- 
designs, the worst-case error with s = 1.5 is calculated using ([Mil and the distance kernel, and 
for spherical t-designs the worst-case errors are calculated by relative results in [12]. From the 
figure we can see that in this case, the computed worst-case errors of approximate spherical t- 
designs and spherical 4.i-designs essentially lie on the same curve, which remains as a conjecture 
that the worst-case errors of both spherical t-designs and spherical 4-designs decay in the same 
speed with respect to the number of points in the case s < 2 on S^. Figure |2(b)| plots the 
worst-case errors for both spherical t-designs and spherical to.i-designs with s = 5.5. From the 
figure we can see that the worst-case errors of spherical fo.i-designs decay faster than the ones 
of spherical t-design with respect to the number of points. 


5 Polynomial approximation on the sphere using spherical tg- 
designs 

5.1 Regularized weighted least squares approximation using spherical t^- 
designs 

In this section we consider the restoration of a continuous function / € (7(8^) from its noisy 
values given at N points Xn = {xi,... ,X 7 v} C by the h — h regularized weighted discrete 
least squares form 


I 

min — 


N L 2k+l 


L 21+1 


j=l 1=0 k=l 1=0 k=l 


(56) 
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where fXj > 0, j = 1,..., N are the weights for each term of the least squares model, A > 0 is 
the regularization parameter, and >0,^ = 0 ,... ,L, k = 1, ... ,2£ + 1 are usually chosen 
with the meaning of certain polynomial operators such as Laplace-Beltrami operator and filtered 
operator OES]. In [23] both a priori choice based physical reason in satellite gravity gradiometry 
problem and a posteriori choice based on reproducing kernel theory are considered to choose 
I3e,k- 

Note that A: = 1,..., 2£ + 1, £ = 0,..., L} is a basis of P^. Problem ([56]) is to find a 

good approximation of / in P/, in the form 


L 21+1 

e=o k=i 

Let the entries of matrix Yl & be 

(YL)i/2+fc = yt,fc(xi), i = l,...,N, £ = 0,...,L, k = l,...,2£ + l, 
and = (/‘^(xi),..., /'^(xat))^. Problem (f56]l can be reformulated as 

min , ;^ll^^(YLa-f‘^)||i + A||Da||i, (57) 

„gR(i+l)2 2 


where 



tiNxN 


Tn 


and D is a diagonal matrix satisfying Op^^/^+k = l^t,k with ^ > 0. For polynomial approxi¬ 
mation on the sphere, an / 2 -regularized weighted least squares model has also been considered 

min ;^||A^(YlQ! - f'^)||iA||Da||2. (58) 

aeR(i+i)2 2 


The regularization of this model is of I 2 norm, which can be seen as a measure of energy. 
It is known that the li regularization has desirable properties in approximation of nonsmooth 
continuous functions. An li regularization term is preferable to be considered here. By choosing 
a suitable penalization term, the I 2 — h regularized model is usually supposed to achieve a 
more sparse solution than the I 2 — h regularized one, which means that the target function is 
approximated by less basis spherical polynomials. Additionally, for functions which are globally 
continuous but locally non-differentiable on the sphere, the I 2 — h regularization is better than 
the I 2 — I 2 regularization. 


Theorem 11. Let X]\f be a spherical t^-design and w be the veetor of weights satisfying 
and ^ with respect to Xjy. Let L > 0. For model [57\) set fij = wj for j = 1,N. Then 


Hl = yIayl = i^L+ir 


(59) 


and has the unique solution 


= max{0, S(,^k - + min{0, S(,^k + >^P(.,k] 

for £ = 0,... ,L, k = l,...,2k + l, where se^k = Ylf=i WiYe,k{^i)f^{^i). 


( 60 ) 
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Proof. Note that when is a spherical t^-design, 

N 

i=l 

= / y£,k{x)Ye',k'{x)duj{x) = 6uOkk', (61) 

JS2 

where the third equality is established by the orthonormality of spherical harmonics. Problem 
(j57l) is strictly convex by the fact that is nonsingular and so it has a unique optimal solution. 
Since A{a) = ^||A 2 (Yia — f '^)||2 is strictly differentiable, by deriving the first optimality 
condition of (j57p and Corollary 1 in [161 Section 2.3], we obtain that its unique optimal solution 
satisfies 

0 G Hlo-Y^W f‘^ + A5(||Da||i), (62) 

where d{-) denotes the subdifferential. By (I6ip which implies and the fact that 

D is diagonal, problem ([621) is separable and thus a is a solution of (1621) if and only if it is a 
solution of 

0 G ~ £ = 0,..., L,/c = 1,..., 2£ + 1. (63) 

Denote = d\a^^k\ and hence —1 < Ti^k < 1- Let be the optimal solution of (f63|) with 
corresponding I and k and hence 


= S£,k - >^l3i,kTl,k with Ti^k G [-1, Ij. 


(64) 


When se^k > ^/3e,k we can set re^k = 1 and obtain 


(^£,k = Si,k — ^I3e,k > 0 , 


which together with fl^^k > 0 satisfies (fHnjl and ill). When s^^k < —^f3£,k similarly we set 
Ti,k = -1 and get 

+ ^Pi^k < 0 , 

which also satisfies (l60]l and (IMl) . Then when S£^k £ [—^Pe,k, ^Pe,k] we set T£^k = J^TZ ^ 
and get that 

“I,A: = 0) 

which also satisfies (|60l) and (1641) . Hence the theorem is proved. □ 

Denote the approximation residual as A{a) = Ylf=iiPL,N{'^j) — /‘^(xj))^. Let a*(A) be 
the optimal solution of (|57D with different regularized parameters A. The following proposition 
indicates that T(q!*(A)) is monotonically increasing with respect to A. 

Proposition 12. Let be a spherical t^-design with t >2L and p.j = Wj for j = 1,..., Y. 
Then A(a*(A)) is increasing in A. 

Proof. Let A, A be given with 0 < A < A and denote the optimal solution of problem (1571) with 
A, A as a*, a* respectively. Denote E{X,a) = A||Da||i and the minimization property of ([57]) 
for A gives 

A{a*) + E{X, a*) < A{a*) + Y(A, a*), (65) 


which implies that 


A{a*) - A{a*) < E{X, a*) - Y(A, a*). 


(66) 
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From (1601) we have 



se,k — ^f3£,k, 
se,k + ^P£,k, 

0 . 


S£,k > ^f3£,k 
S£,k < —^fil,k 
— ^P£,k < 


(67) 


Since > 0, we have \oi*^i.\ = max(0, \s£^k\ — ^P£,k)- Together by the fact that 


L 2£+\ 


E{\a*) = \Y^Y.fit^k\otlkV 


£=0 k=l 


we have 

\a\^k\ = max(0, \s£^k\ - ^/3£,k) < max(0, \si^k\ - A/3^,fc) = |a£,fc|- 
Hence it is obtained that £'(A,d*) < E{X,a*). Together with (f6^ we complete the proof. □ 

5.2 Numerical experiments 

In this snbsection we report the numerical results to test the efficiency of the I 2 — h regularized 
model (f57|) using spherical te-designs. 

Example 5.1. In the first numerical test, the target function is selected as spherical poly¬ 
nomials with degree no higher than L. Obviously using both models (j57p and (1581) the target 
function can be exactly restored when A = 0 and the data is noise free and the optimal values 
of the two models equal to 0 in such case. However, due to the noise in the data vector it is 
necessary to use the regularization models. 

In this experiment we will use the spherical fo.i-designs which is calculated by solving a 
system of nonlinear equation ()12D and the approximate spherical t-designs proposed in m as 
the point set for polynomial approximation. Both the uniform errors and L 2 errors are recorded 
to measure the approximation quality. We choose a large-scaled and well distributed point set 
Xt C to be the test set and use it to estimate the errors. Then the uniform error and L 2 
error of the approximation are estimated by 


\\f - Pl,n\\c{S^) ~ max |/(xi) - pL,Ar(xj)| 

XiGXt 


( 68 ) 


and 



(69) 


where Nt denotes the number of points Xj in Xt- In this experiment, we choose Xt to be an 
equal area partitioning point set |25] with 10^ points. The matrix D in the experiment is always 
selected as for £ = 0,..., L, k = 2i+l, inspired by the Laplace-Beltrami operator, 

see [3]. 

Fig. 5.1 shows the approximation errors using both I 2 — h model ([57]) and I 2 — h model ([58]) 
with different A and different noise scales 6. The noise of the data obeys a uniform distribu¬ 
tion in [—(5, (5]. In this numerical experiment a spherical 37o.i-design with only 514 points, which 
is much less than , is applied to approximate a randomly generated spherical polynomial 

with degree = 18 (The polynomial is generated with all its Fourier coefficients obeying the 
standard normal distribution). The regularization parameter A is chosen from 10“^^ to 
Fig. 5.1 (a)(b) give the errors of the approximation with different A for 5 = 0.1 using the two 
models. From the two sub-figures it can be seen that model (|57p can restore the 18-degree poly¬ 
nomial more accurately than model (|58D . The minimal error with respect to different A can be 
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(a) Uniform Errors with different A for 5 = (b) L 2 errors with different A for 5 = 0.1 

0.1 




S 8 

(c) Minimal uniform errors with different (d) Minimal L 2 errors with different noise 
noise scales S scales S 


Figure 3: Errors for restoring 18-degree polynomial 


achieved at about A = 10 Fig. 5.1 (c)(d) show the errors of the restoration results with differ¬ 
ent noise scales. It can be seen that the model (j57p performs better in each noise scale than (I58p . 

Example 5.2. In the second numerical experiment we test the numerical performance of 
model (1571) using spherical t^-designs and spherical t-designs. We select the Franke function |24j 

/i(x) =f{x,y,z) = 0.75exp(-(9x - 2)V4 - (9y - 2)^4 - {9z - 2)^/4) 

-h 0.75exp(-(9x -h 1)V49 - (9y -h 1)/10 - (9z + 1)/10) 
-h0.5exp(-(9x-7)V4- i9y-3)‘^/4- (9z - 5)V4) 

— 0.2exp(—(9x — Af — {9y — 7)^ — {9z — 5)^), {x,y,z) G 

to be the target function which is not a spherical polynomial but continuously differentiable 
on the whole sphere. We set e = 0.1 and also d = 0.1 in this experiment and the scheme of 
choosing A is the same as in Example 5.1. For spherical to.i-designs we select those point sets 
constructed with possible least points. As is mentioned above, a spherical to.i-design may be 
constructed using less than points. Approximate spherical t-designs proposed in m 

are also applied for comparison. Note that the minimizer of model (1571) has an explicit form 
(jbOp only when t > 2L, so for different t we choose L = [|J. 

For Example 5.2, the approximation errors using both spherical to.i-designs and approxi¬ 
mate spherical t-designs are shown in Eig. 5.2. The X-axis represents the number of points in 
the data sets and the Y-axis represents the minimal uniform errors. From the figure we can see 
that approximation using spherical to.i-designs achieves smaller errors than using approximate 
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(a) Uniform Errors (b) L 2 errors 

Figure 4: Errors for approximating Franke function with different scales of point sets 


spherical t-designs in most cases. Based on the numerical results in Fig. 5.2, the approximation 
quality can be improved with the relaxation of weights using model (|57p . 


Example 5.3. In the third experiment, a continuous but non-differentiable function 


/2 = /l(x) + /cap(x), 

is selected as the target function to approximate, with 

"7rcos“^(xc • x)' 


/cap(x) — 


pcos 


2r 


0 , 


X G C(xc,r), 

otherwise, 


(71) 


(72) 


where p > 0, 0 < r < vr. The function is non-differentiable at the edge of the spherical cap 
C(xc, r). Since the basis functions applied for approximation are spherical harmonic polynomials 
which is globally differentiable on §^, restoration of the edge of C(xc, r) turns to be a challenging 
problem when the data has noise. 

A spherical 37o.i-design with 514 points is used as the data point set in this experiment. 
Other settings in this experiment are S = 0.5, A = 10“^°, 10“^®'®,..., 10^, Xc = (—0.5, —0.5, \/o(5)^ 
r = 0.5 and p = 1- The restorations of /2 using both models ([57)1 and (f58]l are depicted in Fig. 
5.3. Similar with previous experiments, we choose the values of A resulting in minimal uniform 
errors for each model and plot the shape of the restoration function on the sphere. From Fig. 
5.3(c) (d)(e)(f), restoration by model (l57[) is not as smooth as restoration by model (l58[) but has 
smaller errors. And more notably, as highlighted by the rectangle in Fig. 5.3(c)(d), model (|57l) 
restores the non-smooth edges of the spherical cap more accurately than model (I58p . 
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